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摘 要 : 水 化 学 分 类 是 实现 水 资源 可 持续 开发 利用 的 重要 前 提 之 一 。Durov 三 线 图 作为 主要 的 水 化 
学 分 类 图 因 绘 制 复 杂 而 限制 了 其 在 现实 中 的 应 用 , 且 一 些 常见 的 水 化 学 软件 也 不 具备 绘制 Durov 
三 线 图 的 功能 。Durov 三 线 图 包括 原始 型 Durov 三 线 图 及 在 其 基础 上 演变 而 成 的 扩展 型 Durov 三 
线 图 。 通 过 介绍 两 种 Durov 三 线 图 的 差异 、 组 成 特征 和 投影 原理 ,提出 了 在 Durov 三 线 图 中 建立 坐 
标 系 的 方法 ,并 推导 了 将 水 化 学 数据 投影 到 Durov 三 线 图 的 公式 。 在 此 基础 上 ,阐述 了 基于 Matlab 
语言 编程 实现 原始 型 Durov 三 线 图 和 扩展 型 Durov 三 线 图 的 方法 和 流程 ,并 以 公开 发 表 的 数据 绘 


制 了 两 种 Durov 三 线 图 。 基 于 相同 的 数据 比较 发 现 ,相关 水 化 学 软件 绘制 的 原始 型 Durov 三 线 图 
及 一 些 公开 发 表 的 研究 成 果 中 的 扩展 型 Durov 三 线 图 与 本 文 所 用 Matlab 语言 绘制 的 Durov 三 线 图 


相 一 致 ,表明 该 程序 准确 可 行 。 本 程序 既 可 以 实现 原始 型 Durov 三 线 图 的 绘制 ,也 可 以 实现 常见 水 
化 学 软件 不 能 完成 的 扩展 型 Durov 三 线 图 的 绘制 。 该 程序 的 应 用 ,将 有 助 于 Durov 三 线 图 的 高 效 
应 用 ,特别 是 扩展 型 Durov 三 线 图 的 应 用 不 但 可 以 作为 水 化 学 分 类 的 依据 ,同时 还 可 以 指示 某 些 水 


化 学 过 程 。 
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近年 来 ,在 气候 变化 、 人 口 增长 和 经 济 快速 发 展 
的 大 背景 下 ,水 资源 的 总 量 明显 减少 ,然而 人 们 对 有 
限 的 水 资源 的 需求 却 并 没有 因此 减 小 ,水 资源 短缺 
现象 日 益 严重 ,因此 ,对 水 资源 的 可 持续 开发 利用 变 
得 十 分 重要 "0 -9 。 要 实现 水 资源 的 可 持续 开发 利 
用 ,首先 需要 了 解 水 资源 的 水 化 学 类 型 ,进行 水 化 学 
分 类 ,以 便 正确 理解 水 资源 与 环境 相互 作用 的 机 制 ， 
这 对 维持 当地 生态 系统 的 健康 稳定 有 着 极其 重要 的 
MLB, 

水 化 学 分 类 依据 的 是 水 体 中 主要 阴 、 阳 离子 
的 相对 含量 和 比例 关系 。 例 如 前 苏联 舒 卡 列 夫 将 
水 中 六 种 主要 阴 、 阳 离子 中 当量 浓度 大 于 25% 的 
离子 进行 组 合 , 共 分 成 49 类 水 型 "1。 对 于 少量 样 
品 ,采用 该 方法 比较 方便 ,但 是 对 于 大 批量 样品 ， 


文章 编号 : 


为 常见 的 如 Durov 三 线 图 。 该 三 线 图 经 历 了 从 原 
始 型 Durov 三 线 图 (The original Durov Diagram ) 到 
扩展 型 Durov 三 线 图 (The expanded Durov Dia- 
gram ) 的 演变 ,不 但 可 以 用 作 水 化 学 分 类 ,还 可 以 
反映 水 化 学 过 程 ,对 揭示 水 文 地 球 化 学 规律 起 着 
重要 的 作用 ,因而 具有 广阔 的 应 用 前 景 ” 。 然 而 
由 于 Durov 三 线 图 基于 多 个 变量 并 且 投 影 规则 较 
复杂 ,因而 普通 的 水 化 学 软件 都 不 具备 绘制 Durov 
三 线 图 的 功能 ,限制 了 该 图 在 水 化 学 分 类 中 的 应 
用 。 为 了 克服 三 线 图 不 便于 绘制 的 缺点 ， 
CHADHA 利用 基于 直角 坐标 系 的 CHADHA 水 化 
学 图 来 替代 Durov 三 线 图 以 实现 与 Durov 三 线 图 
类 似 的 功能 ” ,国内 也 有 学 者 对 其 进行 了 验 
证 中。 然而 ,虽然 上 述 方法 可 以 避免 Durov 三 线 


该 方法 显得 非常 繁琐 , 且 不 能 直观 的 对 各 样品 的 图 绘制 复杂 的 缺点 ,但 同时 也 失去 了 Durov 三 线 图 
水 化 学 特性 进行 比较 。 水 化 学 分 类 图 克服 了 上 述 可 以 指示 水 化 学 过 程 的 功能 。 此 外 ,Piper 三 线 
缺点 ,可 以 直观 的 表现 不 同 的 水 化 学 类 型 ,其 中 较 图 :具有 与 Durov 三 线 图 类 似 的 水 化 学 分 类 功 
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能 , 且 两 者 指示 的 水 化 学 类 型 本 质 上 是 一 致 的 '” 。 
作为 一 种 常用 的 水 化 学 分 类 图 , Piper 三 线 图 同样 
得 到 了 广泛 的 应 用 ,但 该 图 仍 不 具备 指示 水 
化 学 过 程 的 功能 。 随 着 计算 机 技术 的 发 展 , AL- 
BASSAM 先后 用 BASIC 语言 和 Visual Basic 语言 实 
现 了 Durov 三 线 图 的 计算 机 输出 。 然 而 限于 
操作 系统 和 版 权 等 因素 的 影响 ， 该 程序 无 法 为 科 
研 人 员 所 直接 利用 。 

Matlab 是 由 MathWorks 公司 推出 的 高 性 能 科技 
应 用 软件 , 因 其 可 读 性 强 .调试 简单 等 显著 特点 而 被 
学 术 界 所 熟知 ,尤其 强大 的 图 形 可 视 化 能 力 在 地 学 
研究 中 具有 广泛 的 应 用 前 景 ”: 。 该 软件 不 受 操作 
系统 的 限制 ,可 在 包括 Windows Linux 和 Mac 等 在 
内 的 系统 下 运行 ,方便 了 科研 人 员 之 间 的 交流 。 本 
文 基于 Matlab 语言 编程 实现 了 原始 型 及 扩展 型 Du- 
rov 三 线 图 的 绘制 ,扩展 了 水 化 学 分 类 图 在 实际 科研 
和 工作 中 的 应 用 范围 。 


1 原始 型 与 扩展 型 Durorv 三 线 图 的 
比较 


原始 型 Durov 三 线 图 采用 当量 浓度 计算 百 分 含 
E ,其 由 位 于 中 心 位 置 的 正方 形 投影 区 位 于 正 
方形 投影 区 左 侧 的 阴离子 三 角 图 及 位 于 正方 形 投 影 
区 上 方 的 阳离子 三 角 图 三 部 分 组 成 (图 1) 。 阴 离子 
三 角 图 由 HCO; S047 和 C1 组 成 ,之 所 以 选取 这 三 
种 离子 的 原 Daa E A 的 阴 
离子 ,并 能 反映 水 体 水 化 学 演变 过 程 。 根 据 这 三 种 
Sl Ca r A 组 、 硫 酸 
盐 组 .氧化 物 组 和 混合 组 这 四 种 一 级 类 型 (图 1)。 
次 一 级 的 分 类 由 以 Ca* .Mg 和 Na 这 三 种 阳离子 
产生 的 阳离子 三 角 图 决定 。 为 保证 阴 、 阳 离子 电荷 
平衡 ,Durov 三 线 图 规定 各 阴离子 之 和 为 50% ,各 阳 
离子 之 和 亦 为 50% 。Durov 对 更 早期 的 Durov 三 线 
图 进行 了 改进 1 ,增加 了 pH 和 TDS, 将 其 投影 在 了 
正方 形 的 外 围 (图 1)。 

BURDON 和 MAZLOUM ™! 及 LLOYD "l 进 一 
步 改 进 了 Durov 三 线 图 , 即 扩 展 型 Durov 三 线 图 。 
他 们 将 阳离子 三 角 图 继续 划分 为 三 个 次 一 级 的 单 
元 , 即 两 个 小 三 角形 和 一 个 萎 形 ,然后 将 两 个 小 三 角 
形 从 奏 形 区 域 两 侧 分 别 移动 某 个 单位 的 间隔 , 即 以 
Ca?* 4 25% 和 Na’ 占 25% 线 为 界 , 将 原 阳 离子 三 
角形 一 分 为 三 ;阴离子 三 角 图 也 被 划分 为 两 个 小 三 


角形 和 一 个 次 形 ,然后 将 两 个 小 三 角形 从 菱形 区 域 
两 侧 分 别 移动 某 个 单位 的 间隔 , 即 以 Cl 25% 线 
Al HCO; 占 25% 线 为 界 将 原 阴 离子 三 角形 一 分 为 
三 。 最 后 ,根据 这 些 重新 划分 的 区 域 再 将 原 和 矩形 区 
域 细 化 为 9 个 次 一 级 的 投影 单元 (图 2)。 
Durov 三 线 图 和 原始 型 Durov 三 线 图 相 比 具有 明显 


扩展 型 


的 优势 , 它 提 供 了 一 种 更 好 的 表现 水 化 学 类 型 和 过 


程 的 方式 ,更 能 反映 水 体 的 水 化 学 差异 .相互 关系 和 
离子 来 源 ”” 。 此 外 ,根据 样品 在 9 个 次 一 级 的 投 
影 单元 的 分 布 关系 ,扩展 型 Durov 三 线 图 还 可 指示 
离子 交换 简单 溶解 和 线性 混合 及 道 离子 交换 等 水 


化 学 过 程 9-781 。 
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图 1 原始 型 Durov 三 线 图 


Fig.1 Original Durov Diagram 
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图 2 扩展 型 Durov 三 线 图 及 其 验证 
Fig.2 Expanded Durov Diagram and its validity testing 
注 : 四 :25% <Ca** <50% ;@:0 < Ca** <25% H 0 <Na+ <25%; 
®:25% <Na* <50% ;@:25% < HCO; <50% ;6):0<Cl~ <25% 

and 0< HCO; <25% ;@:25% <CL <50% 


2 Durov 三 线 图 的 绘制 原理 


2.1 原始 型 Durov 三 线 图 的 实现 

(1) 单位 转换 。Durov 三 线 图 以 百 分 含量 表 
示 , 所 以 要 将 各 离子 的 当量 浓度 值 转换 为 对 应 的 百 
分 含量 。 需 要 注意 的 是 , 阴 、 阳 离子 各 占 50% , 而 非 
100% 。 

(2) 建立 坐标 系 , 绘 制 边框 。 首 先 , 定 义 阴 、 阳 
离子 三 角 图 交点 位 置 为 坐标 原点 (图 3) ,并 规定 三 
角形 边 长 为 50, 则 可 求 得 阴 、 阳 离子 三 角 图 及 和 矩形 
区 域 各 顶点 的 坐标 ;其 次 ,根据 求 得 的 坐标 点 ,用 
plot 命令 绘制 Durov 三 线 图 各 部 分 (阴阳 离子 三 角 
图 及 和 矩形 区 域 ) 的 边框 。 

(3) 坐标 转换 及 数据 投影 。 在 已 知 某 样品 各 高 
子 百 分 含量 的 情况 下 ,将 其 直接 投影 到 三 线 图 坐标 
系 中 并 不 容易 ,简单 的 做 法 是 先 将 各 离子 百 分 含 量 
转换 为 直角 坐标 系 下 点 的 坐标 ,然后 直接 投影 到 直 
角 坐 标 系 中 。 

假定 样品 为 A, 各 离子 百 分 含 量 已 知 , 设 其 在 阳 
离子 三 角 图 中 的 坐标 为 xm ,yi,A 的 Ca Boat 
Ha, Na! ESERE b, M Mg 的 百 分 含 量 ; 
50% -ww -bi; 设 A 在 阴离子 三 角 图 中 的 坐标 为 x, 
yy, 其 Cl 百 分 含 量 为 a,,HCO; BASEEN b, M 
S01 的 百 分 含 量 为 50% - a, -5,; 设 A TERETE Sak 
的 坐标 为 xs ,” 。 根 据 图 3 阴影 部 分 的 三 角 了 水 数 关 
系 , 经 简化 后 可 以 得 出 样品 A 在 直角 坐标 系 下 投影 
到 阳离子 三 角 图 中 对 应 的 坐标 值 : 

x, =0.5 x (50% -al +b,) (1) 

y, =0.5% xtg 60° x (50% -a,-b,) (2) 

同 理 ,样品 A 在 直角 坐标 系 下 投影 到 阴离子 三 角 图 
中 对 应 的 坐标 值 为 : 

x, = -0.5 xtg 60° x (50% -a,-b,) (3) 

y, = -0.5 x (50% -a, +b,) (4) 

则 样品 A EHARA TR BI) HBT DX Sak POY Ni 

的 坐标 值 为 : 


%3 =X] (5) 


y3 = Y2 (6) 
最 后 ,将 a,b, ay b, 各 值 分 别 带 入 到 式 (1)、(2)、 
(3) (4) (5) 和 (6) 中 求 出 样品 投影 到 原始 型 Du- 
rov 三 线 图 对 应 区 域 的 坐标 值 ,并 用 plot 命令 将 其 绘 
制 在 相应 位 置 。 


图 3 原始 型 Durov 三 线 图 的 坐标 转换 


Fig.3 Coordinate transformation processes 


for the original Durov Diagram 


2.2 ”扩展 型 Durov 三 线 图 的 实现 

扩展 型 Durov 三 线 图 的 绘制 方法 与 原始 型 Du- 
rov 三 线 图 大 体 相 同 ,不 同 点 在 于 扩展 型 Durov 三 线 
图 将 阴 、 阳 离子 三 线 图 进一步 细 分 ,相应 的 矩形 投影 
区 也 被 细 化 。 具 体 说 来 ,在 阳离子 投影 区 ,根据 阳 离 
子 百 分 含量 的 分 布 范围 ,其 被 投影 到 三 个 不 同 的 区 
域 。 当 25% < Ca <50% 时 ,样品 将 被 投影 到 位 于 
阳离子 投影 区 左 侧 位 置 的 一 个 三 角形 区 域内 ( 即 图 
2 RAND) ; 当 0 < Ca <25% H 0 <Na* <25% Hf, 
样品 将 被 投影 到 位 于 阳离子 投影 区 中 间 位 置 的 净 形 
区 域内 ( 即 图 2 所 示 @); 当 25% <Na* <50% R}, 
样品 将 被 投影 到 位 于 阳离子 投影 区 右 侧 位 置 的 三 角 
形 区 域内 ( 即 图 2 所 示 @)。 扩 展 型 Durov 三 线 图 的 
实质 是 将 原始 型 Durov 三 线 图 的 阳离子 三 角 图 切割 
为 三 部 分 并 进行 平移 ,使 其 外 形 类 似 于 Piper 三 线 
图 。 同 理 ,阴离子 投影 区 也 被 切割 平移 为 三 部 
分 , 即 当 25% < HCO; <50% 时 ,样品 将 被 投影 到 位 
于 阴离子 投影 区 上 部 的 三 角形 区 域 ( 即 图 2 所 示 
@);*40< Cl <25% HO< HCO; <25% 时 ,样品 
FR BN) (iF BA PS FE ae DB EY BT DX ak (BI 
图 2 MRO); 4 25% < Cl <50% 时 ,样品 将 被 投 
影 到 位 于 阴离子 投影 区 下 部 的 三 角形 区 域 ( 即 图 2 
ERO) 。 因 此 ,在 将 样品 投影 到 扩展 型 Durov = 
图 的 过 程 中 ,除了 将 样品 按照 原始 型 Durov 三 线 图 
的 方法 求 坐标 ,还 要 将 求 得 的 坐标 值 进 行 坐标 平移 。 

设 样 品 为 A, 且 其 基本 参数 仍然 使 用 原始 型 
Durov 三 线 图 中 所 用 的 方式 表示 ,坐标 平移 距离 这 
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里 取 12.5, 则 样品 A 在 直角 坐标 系 下 投影 到 阳离子 


三 角 图 中 对 应 的 坐标 值 为 : 
当 25% <Ca* <50% 时 ， 
x, =0.5 x (50% -ai +b,) (7) 
yı =0.5 x tg 60° x (50% -a, -b,) (8) 
当 0 < Ca <25 H 0 <Na* <25 H}, 
x, =0.5 x (50% -a, +b,) +12.5 (9) 


yı =0.5 xtg 60° x (50% -a,-6,) (10) 

x4 25% <Na <50% 时 ， 
x, =0.5 x (50% -a +b,) +25 (11) 
yı =0.5 xtg 60° x (50% -a,-b,) (12) 


样品 A 在 直角 坐标 系 下 投影 到 阴离子 三 角 图 
中 对 应 的 坐标 值 为 : 
当 25% <Cl <50% 时 ， 


x, = -0.5 xtg 60° x (50% -a,-b,) (13) 
yı = -0.5 x (50% -a, +b,) -25 (14) 

“O<Cl <25% H 0< HCO; <25% 时 ， 
x, = -0.5 xtg 60° x (50% -a,-b,) (15) 
y = -0.5 x (50% -a, +b,) -12.5 (16) 

™4 25% < HCO; <50% 时 ， 

x, = -0.5 xtg 60° x (50% -a,-b,) (17) 
y, = -0.5 x (50% -a, +b,) (18) 
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3 Durov 三 线 图 的 编程 流程 及 其 验证 


3.1 编程 流程 及 误差 分 析 

编程 实现 Durov 三 线 图 输出 的 关键 包括 求 坐标 
和 绘 边框 两 部 分 (图 4) 。 样 品 的 阴 、 阳 离子 含量 通 
常用 mg + L 表示 , 求 坐标 是 将 该 单位 的 数据 首先 
转换 为 当量 浓度 (meq* L”) ,再 分 别 在 阴 、 阳 离子 
总 浓度 为 50% 的 条 件 下 , 求 出 对 应 各 离子 的 百 分 含 
量 , 然 后 根据 本 文 第 三 部 分 给 出 的 公式 计算 样品 在 
直角 坐标 系 下 的 对 应 坐标 。 绘 边框 是 在 确定 坐标 原 
点 后 ,用 plot 命令 绘制 出 组 成 Durov =28 A AY BA 、 阳 
离子 三 角形 投影 区 域 . 稳 形 投影 区 域 和 算 形 投影 区 
域 等 的 过 程 。 最 后 ,用 plot 命 令 将 带 有 坐标 的 样品 


数据 输入 
单位 /mg + L') 
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原点 
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图 4 Durov 三 线 图 编程 流程 


Fig.4 Program flow of plotting Durov Diagram 


点 投影 到 Durov 三 线 图 的 相应 位 置 。 

在 整个 运算 过 程 中 误差 的 出 现 主要 体现 在 以 下 
几 个 方面 :首先 ,在 将 离子 浓度 从 mg ' 工 ”转换 为 当 
量 浓度 时 ,相对 原子 质量 有 效 数字 位 数 的 选取 会 影 
响 到 结果 的 精度 ,本 文采 用 小 数 点 后 保留 3 位 有 效 
数字 的 原则 ;其 次 ,在 整个 转换 过 程 中 的 数据 运算 及 
其 后 的 结果 所 采用 的 数值 类 型 也 会 影响 到 最 终 绘 图 
的 精度 ,本 文采 用 Matlab 中 精度 最 高 的 双 精 度 浮 点 
数 ( 即 Matlab 中 的 double 数值 型 ) ;再 次 ,将 离子 的 
当量 浓度 转换 为 百 分 含 量 及 计算 样品 在 直角 坐标 系 
下 的 坐标 时 数值 类 型 的 选择 ,同样 会 影响 绘图 的 精 
度 , 本 处 仍 采 用 双 精 度 浮 点 数 。 综 上 所 述 , 本 文 所 采 
用 的 有 效 数字 位 数 和 数值 类 型 可 保证 绘图 结果 具有 
较 高 的 精度 。 
3.2 与 相关 研究 成 果 的 比较 

为 了 验证 程序 的 正确 性 ,采用 经 典 文献 LLOYD 
和 HEATHCOTE |*! 48 6 章 第 四 节 Durov 三 线 图 相关 
部 分 的 数据 绘图 ,数据 见 表 1 ,所 绘制 的 原始 型 
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图 5 对 原始 型 Durov 三 线 图 的 验证 


Fig.5 Validity testing for the original Durov Diagram 


#1 LLOYD 和 HEATHCOTE "! 文中 阴 、 
阳离子 数据 / meq - L“ 
Tab.1 Data of anions and cations in LLOYD and 
HEATHCOTE |*! / meq - L`! 


Ca2+ Mg?* Nat K* HCO; S04- Cl- pH TDS 


q 2a 


1 7.58 0.48 0.74 0.07 6.61 0.60 1.30 2.70 660 
2 0.40 1.28 8.09 0.18 5.34 0.77 3.47 10.40 710 
3 4.24 2.88 11.74 0.15 7.84 0.60 10.89 3.40 1 190 


Durov 三 线 图 (图 5) ,扩展 型 三 线 图 (图 2)。 由 于 相 
关 文 献 [8] 较 早 ,其 原始 型 Durov 三 线 图 不 同 于 当前 
主流 的 阴 、 阳 离子 分 布 格局 , 即 其 阴离子 投影 区 位 于 
正方 形 投影 区 的 上 方 而 阳离子 投影 区 位 于 正方 形 投 
影 区 的 左 方 , 这 与 当前 主流 的 阴阳 离子 分 布 格局 相 
反 。 通 过 对 比 LLOYD 和 HEATHCOTEI5 原始 型 Du- 
rov 三 线 图 与 图 5 数据 的 投影 位 置 ,发 现 两 者 本 质 上 
是 一 致 的 。 此 外 ,LLOYD 和 HEATHCOTE 号 文中 
的 扩展 型 Durov 三 线 图 与 当前 主流 模式 一 致 , 即 阳 
离子 投影 区 位 于 正方 形 投影 区 上 方 、 阴 离子 投影 区 
位 于 正方 形 投影 区 左 方 ,经 与 图 2 对 比 发 现 两 者 表 
示 的 样品 位 置 完全 一 致 。 以 上 结果 表明 本 程序 中 的 
原始 型 Durov 三 线 图 和 扩展 型 Durov 三 线 图 均 正确 
可 用 。 
3.3 与 其 他 水 化 学 分 析 软 件 的 比较 

目前 可 绘制 Durov 三 线 图 的 商业 软件 并 不 多 ， 
且 常 见 的 如 AquaChem 和 AqQA 等 仅 能 绘制 原始 型 
的 Durov 三 线 图 而 无 法 绘制 扩展 型 Durov 三 线 图 。 
为 了 验证 本 文 方法 的 可 靠 性 ,仍然 使 用 LLOYD 和 
HEATHCOTE 的 数据 ,并 采用 2014. 2 版 本 的 
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AquaChem 软件 绘制 原始 型 Durov 三 线 图 ,结果 见 图 
6。 通 过 比较 发 现 用 AquaChem 软件 绘制 的 原始 型 
Durov 三 线 图 (图 6) 和 本 文 用 Matlab 编程 绘制 的 原 
始 型 Durov 三 线 图 (图 5) 并 无 本 质 区 别 ,主要 差异 
在 于 pH 部 分 。AquaChem 软件 中 pH 的 刻度 值 自 上 
而 下 被 设置 为 递减 ,而 本 文 用 Matlab 绘制 的 pH 部 
分 为 保持 和 经 典 文献 相 一 致 而 采用 刻度 值 自 上 而 下 
递增 的 设置 。 商 业 软 件 的 优势 在 于 操作 简单 ,但 因 
其 代码 未 开源 ,因而 不 能 进行 二 次 开发 。 相 反 ,使 用 
Matlab 编程 实现 的 Durov 三 线 图 可 根据 需求 进行 相 
应 的 更 改 , 例 如 对 原始 型 Durov 三 线 图 代码 做 相应 
编辑 后 即 可 得 到 扩展 型 Durov 三 线 图 ,这 是 商用 软 
件 无 法 实现 的 。 
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图 6 AquaChem( 版 本 2014. 2) 绘 制 的 原始 型 
Durov 三 线 图 
Fig.6 Original Durov Diagram ploted by AquaChem 
( Version 2014.2) 
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Durov 三 线 图 作为 重要 的 水 化 学 分 类 图 ,长 久 
以 来 都 因为 绘制 复杂 而 使 其 实际 应 用 受到 了 限 
制 ,特别 是 扩展 型 Durov 三 线 图 无 法 用 常见 的 水 化 
学 软件 绘制 。 本 文 在 详细 介绍 原始 型 及 扩展 型 
Durov 三 线 图 的 组 成 和 绘制 原理 的 基础 上 ,以 当前 
科研 工作 者 常用 的 Matlab 语言 为 工具 ,实现 了 该 
图 的 计算 机 编程 绘制 ,并 将 其 与 常见 商用 软件 给 
制 的 Durov 三 线 图 及 公开 发 表 的 研究 成 果 中 的 
Durov 三 线 图 进行 了 比较 ,表明 该 程序 准确 可 行 。 
该 程序 不 但 可 用 实现 相对 简单 的 原始 型 Durov 三 
线 图 的 输出 ,同时 还 可 以 实现 一 般 商 业 软 件 难 以 
绘制 的 扩展 型 Durov 三 线 图 。 该 方法 的 应 用 将 有 
助 于 水 化 学 分 类 的 高 效 实现 ,同时 也 可 以 作为 指 
示 某 些 水 化 学 过 程 的 工具 。 
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Plotting Durov Diagram based on Matlab 


REN Xiao-zong', LIU Min', ZHANG Ying-zhen', HE Zhi-ming', ZHU Bing-qi? 
(1 Taiyuan Normal University ,School of Geographical Sciences , Jinzhong 030619 , Shanxi , China ; 
2 Key Laboratory of Water Cycle and Related Earth Surface Processes , Institute of Geographic Sciences and 
Natural Resources Research , Chinese Academy of Sciences , Beijing 100101 , China) 


Abstract; Hydrochemical classification of water resource is one of the vital preconditions for water resource sus- 
tainable utilization. There are numerous methods with regard to hydrochemical classification for water resource , but 
the most effective and simplest method is the diagram method. As one of the most important diagram methods for 
hydrochemical classification ,the Durov diagrams (including both the original type of Durov Diagram and the ex- 
panded Durov Diagram) are used widely not only because it can reveal the hydrogeochemical characteristics of wa- 
ter, but also because it can show the hydrochemical processes occurring within different hydrogeological systems. 
However , comparing with other diagrams (such as the Piper Diagram) ,these diagrams are very complex and diffi- 
cult to plot. Few commercial software packages are designed to plot Durov diagrams, which limits researchers to use 
the diagrams. In this paper, we introduced both the original type of Durov Diagram and the expanded Durov Dia- 
gram ,and we also compared the similarities and differences in them. To better understand the diagram, we interpre- 
ted the composition of the original type of Durov Diagram, especially the part of cation triangle , anion triangle and 
the square field. We also interpreted the composition of the expanded Durov Diagram, in which each of the cation 
and anion triangles is split into two subtriangles and a diamond, while the square field is split into nine fields. We 
stressed how to establish a coordinate to plot the original type of Durov Diagram and the expanded Durov Diagram 
respectively. And we deduced formulas to transform the hydrochemical data so that they can be projected into both 
the original type of Durov Diagram and the expanded Durov Diagram. Using these principles, we wrote programs for 
the output of both the original type of Durov Diagram and the expanded Durov Diagram based on Matlab. To test the 
validity of the programs ,we plotted the original type of Durov Diagram with published data based on our Matlab pro- 
grams. And we also plot the original type of Durov Diagram with the same data based on commercial software. The 
results showed that our program was as correct as the commercial software , which indicated our program is feasible. 
However , because none of commercial software packages we found has the function of plotting the expanded Durov 
Diagram ,we compared the diagram plotted by our Maltlab program with diagram appearing in some classical litera- 
ture that used the same data. The results also confirmed that our Matlab program is reliable. Both the original type of 
Durov Diagram and the expanded Durov Diagram plotted by Matlab will help researchers to interpret hydrochemical 
data in graphical methods conveniently. Especially, the expanded Durov Diagram provides a method to indicate 
hydrochemical processes. 


Key words; hydrochemistry; the Durov diagrams; Matlab; plotting 


